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Abstract 



We investigate an event excess in the reaction e'^p — > e^X with large x and 
. high observed at HERA and show that the events could naturally be in- 

r> I terpreted as a signature of the production of the scalar top quarks (stops) 

' in a supersymmetric model with i?-parity breaking interactions. The HERA 

events are characterized by the broad mass distribution and in fact it can be 
simulated by our specific scenario if we consider almost degenerate two mass 
eigenstates ti and t2 of the stops. 



1 



Recently both the HI and ZEUS Collaboration reported that an event excess in the 
deep inelastic scattering (DIS) e~^p — > e'^X with large x and high had been observed 
H, H at HERA. The total data sample analyzed at the HI and ZEUS corresponds to 
respectively accumulated luminosity of 14.2pb^^ and 20.1pb~^ in positron (27.5GeV) - 
proton (820GeV) collisions. For high (> IS.OOOGeV^) region 12 events are observed 
at both the HI and ZEUS, while the Standard Model expectation are 4.7 it 0.8 and 8.7 it 
0.7, respectively. Combining these results we get Nobs = 24 {(Jobs ~ 0.71pb) against N^xp 
= 13 ± 1 {aexp ^ 0.38pb). 

The event excess is characterized by not only a high four momemtum transfer 
but also a large Bjorken parameter x ^ 0.3. Interestingly, the events seem to cluster at 
invariant mass M = 200 ~ 220GeV, where M is directly determined from x hy M = y/xs. 
However, while the excess events of the HI seem to cluster at M ~ 200GeV, those of the 
ZEUS have a broader distribution around M ~ 220GeV. Statistical consistency between 
the HI and ZEUS data was discussed in ref.j^, where a possible single resonance structure 
was assumed. Here we will simply combine both data since each statistical significance is 
still not so high enough. 

Some possible interpretations of the events have already been given ; 
(A) four fermion contact interaction §, 11| 



(B) leptoquark production ^, 0, 11, 

(C) scalar quark production in the supersymmetric (SUSY) model with R-parity breaking 

(RB) interactions g |, |, |, |lO|, 0, H 



In this letter we investigate one of scenarios of the category (C), the scalar top (stop) ||T_ 
production in the RB SUSY model ||^, |5|, |^, 0. For we have first proposed the scenario at 
HERA experiments long since |14, |l^]. On the other hand, possibilities of the search for 
the squarks in the first or second generation at HERA have originaly been discussed in 
ref.H, PI. 



The discussion is based on the minimal SUSY standard model (MSSM) with an RB 

interaction 

L = X[j^.{u]ldkPLe - dkB^PLUj) + h.c, (1) 

where Pl,r read left and right handed chiral projection operators, respectively. The 
interaction Lagrangian (|l[) has been derived from the general RB superpotential |]l^]; 

Wjf = XijkLiLjE^ + X'ijkUQjD^k + X'ljkU'iD^jD%, (2) 

where i,j,k are generation indices. The first two terms violate the lepton number L 
and the last term violates the baryon number B. If we want to explain such unresolved 
problems as (i) the cosmic baryon number violation, (ii) the origin of the masses and the 
magnetic moments of neutrinos and (iii) some interesting rare processes in terms of the 
L and/or B violation, i?-parity breaking terms must be incorporated in the MSSM. The 
coupling Eq. (|l|) will be most suitable for the ep collider experiments at HERA because 

the squark UjL or d^R will be produced in the s-channel in e~^-q sub-processes. 

+ dk ^ uJl ^ + dk (3) 
+Uj ^ dkR + Uj. (4) 
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The upper bounds on the couphng constants have already been settled by some 



experiments, neutrinoless double beta decay [^], charged current universality [1^, 2C] 



atomic parity violation (APV) |18| , pOj , 21| and Vf. mass [22|. Among nine candidate 
coupling constants A'^^^ some are severely constrained by these experiments. To explain 
the observed cross section (~ 0.7pb) as well as to circumvent the experimental constraints, 
we should consider the two kinds of processes (3) in which a valence quark in the proton 
participates, 

e^ + d^cl^e^ + d (5) 
e+ + d ^ tl ^ e+ + d. (6) 

In fact, the stringest upper bounds on the corresponding coupling constants coming from 
the APV experiments are not so severe, 

A'i2i,A;3i<0.13 (7) 

for fn~^, ~ 200GeV. Note that the scalar charm (scharm) cl and the stop cannot 
couple to any neutrinos via i?-breaking interactions. That is, no event excesses are ex- 
pected in the charged current process ep — > uqX in the squark scenarios. This is a unique 
property of the squark scenarios which could be useful for us to distinguish the squarks 
from leptoquarks. 

Here we pay attention to a fact that the stops {t^, t^) are naturally mixed each other 
|p!3| due to a large top quark mass [^] and the mass eigenstates (ti, t2) are parametrized 
by a mixing angle 8t, 

ti = ticosOt — t2sm6t. (8) 
In this case the interaction Lagrangian (||) for the stops is written by 

L = A'lgi (cos OttidPLe - sin ett^APLe) + h.c. (9) 

It should be emphasized that only one scharm cl could be coupled to e^d in the scenario 
with A'^21 7^ 0- 

Before discussing numerical results, we examine the decay modes of the stop. In the 
MSSM, the stop lighter than the other squarks in the 1st and 2nd generations and the 
gluino can decay into various final states : 

t - 





(a) 


bW, 


(b) 


Wb 


(c) 


b£u 


(d) 


bvl 


(e) 


bWZk 


(f) 


bfJZk 


(g) 


c Z\ 


(h) 


e d. 


(i) 



where Z^ {k = \ ^ 4), Wi(i = 1,2), v and I, respectively, denote the neutralino, the 
chargino, the sneutrino and the charged slepton. (a) ~ (h) are the i?-parity conserving 
decay modes, while (i) is only realized through the RB couplings (|9|). 
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If we consider the RB couphng with X[^i > 0.01, the decay modes (d) to (h) are 
neghgible due to their large power of a arising from multiparticle final state or one loop 
contribution. Moreover, in the present case (m^^ ~ 200GeV) the mode (a) will be kine- 
matically suppressed. Then only two body decay modes (b), (c) and (i) are left for our 
purpose. 

In Fig.l we show the X'l^^ dependence of the total cross section for the process, 

e+p ^ hX e+gX, (10) 

where we take m~ = 210GeV, 9t = and > 15,000GeV^. The horizontal lines corre- 
spond to the observed cross section at HERA, aobs = (0.71 it 0.15) pb. We find the HERA 
data could be explained by A'^gi ~ 0.04, 0.06 and 0.13 for Br(ii e+d) = 1.0, 0.5 and 0.1, 
respectively. This is almost the same result with that given by Altarelli et al. [^. From 
the Fig.l we see that 0.03 ~ X'l^i < 0.13 and Br(ti — > e~^d) ~ 0.05 are needed for us to 
obtain the required amount of the cross section of about 0.7pb. These facts immediately 
lead us to another important consequence, Ttot (t) ~ 1.5GeV, because 

r(fi^ e+d) = ^^H^^mr ^0.07 GeV (11) 

and Te^tra < (1 - 0.05)/0.05 T{ti ^ e+d), where Textra = Ttot(t) - T{ti ^ e+d). Note 
that results obtained here are not altered even if we consider the stop mixing 9t ^ 0. 
Moreover, they are applicable to the scharm scenario when we simply replace {ti — > cl, 

-^131 ~^ -^'121 )• 

Br(t e+d) = 1.0 can easily be realized for m~ < mf, + and m~ < mw + 

g. In this case it is difficult to discriminate the stop from one of the leptoquark S1/2 with 
the charge Q = |. Such a stop (or leptoquark) could be searched at Tevatron through the 
production mechanism, 

pp^tiX ^ eeX (12) 

Recently -DO group has reported preliminary bounds on a leptoquark mass as Mlq > 
194GeV for Br(LQ e+q) = 1.0 [|4|. As has been pointed out by Altarelh et al. [|] 
and Dreiner et al. |^], detailed analyses at DO and a joint analysis of the DO and CDF 
dielectron data will be desirable for discovering or rejecting the stop (or leptoquark) with 
Br(t ^ e+d) = 1.0. 

On the other hand, if the Wi or b is light, i.e., ~ ~ "^VKj 

Br(t bWi) or Br(t — > bW) can compete with Br(t e+d). In this case, we should take 
into account of the processes 

ep bWiX, (13) 

ep bWX (14) 

It is expected that the detectable cross sections a ~ 0.1 pb for heavy stop with mass 
mr~ 250 GeV for e+ beams [15|. In our model the LSP, the lightest neutralino Zi, and b 



We cannot easily set Br(cL ^ d) — 1.0 because cl can decay into cZi via R-conserving interaction 
with a large branching fraction for a wide range of SUSY parameter space Q. 
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possibly decay into i?-even particles via only non-zero RB coupling X'l^i- A typical decay 
chains will be 

ep bWiX b{£uZi)X b{iu{bdu))X, (15) 

ep^bWX {bZi)WX ^ {b{bdi^)){eu)X (16) 

^ {ued)WX {ved){lv)X (17) 

Thus, a possible typical signature of the stop production ep — > bWiX or ep bWX 
would be 26-jets+jet+lepton+ or jet+lepton+ ^t- One of the signals to be detected 
at HERA is characterized by the high Pt spectrum of muons. When such event excess 
could be observed, the stop could clearly be distinguished from the leptoquark Si/2- 

Next we investigate the and M distributions of the expected events. In this analysis, 
we use the and Mg determined from energies and scattering angles of the final positron 
in combining the HI and ZEUS data which are explicitly presented for the selected events 
refs.|,|. 

Qg and Mg distributions of the expected number of events together with the experi- 
mental data are shown in Figs. 2 and 3. Figure 2 corresponds to the case that only the 
lighter stop ti with mass mr = 210GeV is produced and the heavier one t2 is decoupled 
{9t = 0). While the event excess at high Qe is reproduced, the broad Me distribution 
cannot be well explained. It should be noted that the total decay width of the single stop 
is smaller than about l.SGeV as mentioned in the event rate analysis. Consequently, the 
observed broad distribution would not be well explained by the single stop resonance 
scenario or by equivalent scharm cl scenario. 

In Fig. 3, on the other hand, we consider the case that two stops are almost degenarate 
in mass {mj^ , m-- ) = (205GeV, 225GeV) with the finite mixing angle 9t = 0.95. In this case 
double resonance peaks could overlap each other and simulate the broad distribution. 
The event excess at high of the data is also reproduced by this scenario. Difference 
between the one stop and the two stops scenario would be found in the very high (> 
40,000GeV^) region. For simplicity, we take Br(ti^2 e'^d) = 1.0 in the calculation of 
Figs. 2 and 3. It should be noted that the broader M distribution can be obtained for 
Br(ti,2 e+d) < 1.0 than for Br(?i,2 ^ e+d) = 1.0. 

It has been pointed out that the ALEPH 4-jet events [^] are possibly explained by the 
production mechanism e^e~ — > 6^6^ cdcd in the scharm scenario with light sleptons 
mj_ ~ 60GeV [26|. We have discussed the 4-jet events in the stop scenario [27|. In this 
case possible mechanism is e+e~ — > i^i^g — > bdbd. Recent analysis of the 4-jet events 
seems to reveal no 6-quark enhancements and pair production of equal mass particles will 
be disfavoured. If ths is the case, our stop scenario may have no relation to the ALEPH 
4-jet events. 

We have investigated a possible scenario to explain an event excess in the reaction 
e~^p e~^X with large x and high observed by the HI and ZEUS Collabolation at 
HERA by the resonance production of the stops with an i?-parity breaking interaction in 
the framework of the MSSM. The HERA events are characterized by not only high but 
also the broad M distribution. It can be simulated by our specific scenario if we consider 
almost degenerate two mass eigenstates ti and t2 of the stops. 
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Our scenario would be confirmed or rejected at HERA or Tevatron through the search 
for the high Pt muon events (HERA) or the high Pt ee events (Tevatron). Moreover, in 
our favourable scenario, mj^ ~ 210GeV, one of sbottom 6^ could be light as 120GeV ~ 
150GeV and it can decay into v^d via the RB interactions proportional to It will be a 
good target of LEP2 or future linear colliders. Detailed studies of the sbottom production 
at e^e" colliders in our specific scenario will be reported soon. Needless to say, for our 
purpose it would be highly desirable to carry out the high luminosity run at HERA for 
making sure the broad M distribution. 
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Figure Captions 



Figure 1: X'l^i dependence of the total cross section a(e'^p tiX e'^qX). We take 
m^^ = 210GeV, 9t = and > 15,000GeV^. The horizontal lines correspond to 
the experimental value at HERA, aobs = (0.71 ± 0.15) pb. 

Figure 2: and Mg distribution of the expected number of events together with the 
experimental data. We take mj^ = 210GeV, 9t = 0, X'l^i = 0.04 and integrated 
luminosity L = 34pb~^. For Me distribution > 15,000GeV^ and ye > 0.4 are 
adopted. Dashed line corresponds to the SM expectation. 

Figure 3: Qg and Mg distribution of the expected number of events together with the 

experimental data. We take m~ = 205GeV, mj^ = 225GeV, Ot = 0.95, \[^^ = 0.045 
and integrated luminosity L = 34pb^^. For Mg distribution > 15,000GeV^ and 
ye > 0.4 are adopted. Dashed line corresponds to the SM expectation. 



9 




Fig.1 




c 
w 

> 
LU 




0\ (Gey ) 



200 220 240 

(GeV) 



260 



Fig.2 





200 220 240 

(GeV) 



260 



